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حکیده 


خشکسالی یکی از پدیده‌های طبیعی در اقلیم‌های مختلف با ویژگی تکرار شونده است که به علت داشتن 
ویژگی‌هایی نظیر شروع خاتمه شدت و فراوانی» دارای پیچیدگی‌هایی در مقایسه با سایر مخاطرات 
طبیعی است و گستره وسیعی از ابعاد محیط طبیعی و جوامع انسانی را تحت تأثیر خود قرار می‌دهد؛ 
ازجمله بخش‌های اقتصادی, اجتماعی و زیست‌محیطی. در این پژوهش سعی شده تا مدلی برای پیش‌بینی 
ریسک خشکسالی. ارائه و اقدامات ازم به‌منظور کاهش اثرات ناشی از آن انجام شود. بدین منظور از 
شاحص‌های خشکسالی هواشناسی ٩۳۳]‏ و 618191 و شاحص خشکسالی هیدرولوژیک (5111) استفاده شد. 
به‌منظور بررسی تأثیر تغییر اقلیم بر این پدیده دو سناریوی انتشار خوش‌بینانه و بدبینانه مورد بررسی قرار 
گرفت. دوره پایه این مطالعه ۱۹۸۲-۲۰۱۵ و دوره پیش‌بینی ریسک خشکسالی سال‌های ۲۰۲۰-۲۱۰۰ 
است. دو روش مدل‌سازی 80۴ و 627690100 00۵272۳0۵۲16 10699 برای مدل‌سازی ریسک 
خشکسالی استفاده شد. نتایج بدست‌آمده از اين مطالعه حاکی از کاهش نسبی در وقوع دفعات خشکسالی 
در ۲۷ ساله آینده نزدیک در زیر حوضه آبریز افین است. همچنین مقایسه روش‌های مختلف نشان از 
مناسب بودن مدل رگرسیون ناپارامتری برای پیش‌بینی ریسک خشکسالی آینده را داشت. 

کلیدواژه‌ها: مخاطره خشکسالی. آسیب‌پذیری خشکسالی. مدل‌های اقلیمی»136۳4.5 10۳8.5 
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۱۱/۸ جغرافیا و مخاطرات محیطی شمارة سی و سوم 


۱- مقدمه 

خشکسالی یکی از پدیده‌های طبیعی و مکرر اقلیمی در زمین است که اثرات متضوعی بر سایر بخش‌های 
اقتصادی, اجتماعی و زیست‌محیطی ایجاد کرده است. پاندی و همکاران" در سال ۲۰۱۲ در مطالعه‌ی خحود اعلام 
کردند که خشکسالی بر اساس تأثیر مستقیمی که بر زندگی انسان‌ها دارد بالاترین رتبه در میان سایر مخاطرات طبیعی 
را دارد. تأثیر خشکسالی بر کشاورزی, دسترسی به آب و فعالیت‌های اجتماعی و اقتصادی در هر منطقه به عوامل 
فیزیولوژیکی, محیطی و اجتماعی آن منطقه بستگی دارد؛بنابراین آسیب‌پذیری خشکسالی بستگی به عوامل متعددی 
نظیر ویژگی‌های توپوگرافی» خاک زمین مورد استفاده. توسعه منابع آب. بهره‌برداری از آب‌های زیرزمینی و نیاز آب 
منطقه‌ای برای فعالیت‌های داخلی» صنعتی و کشاورزی دارد (آرنل /1۹۹۹). ارزیابی تأثیرات خحشکسالی در بخش‌های 
مختلف دشوار است» زیرا تأثیرات می‌تواند به صورت منطقه‌ای یا محلی باشد. همچنین» خشکسالی‌ها ممکن است 
برای مدت طولانی (یک یا چند سال) و یا مدت کوتاهی (چند هفته) اتفاق بیفتد (وو و وایت "۲۰۰8). حشکسالی 
یکی از بزرگ‌ترین خحطرات طبیعی است که تأثیرات منفی گسترده‌ای بر محیط اقتصاد و جامعه دارد (بچمایر و 
همکاران" ۲۰۱۷). حشکسالی به دلیل کاهش بارندگی در یک دوره مشخص و کارآمد. منجر به کمبود آب می‌شود و 
خشکسالی هیدرولوژیکی به عنوان کاهش آب موجود در تمام اشکال آن تعریف می‌شود (ما و همکاران» ۲۰۱۵). 
پراساس تعریف برنامه استراتژی‌های بین‌المللی برای کاهش بلایا (۲۰۰۷» ریسک خشکسالی حاصل قرار گرفتن یک 
منطقه در برابر یک واقعه خشکسالی با شدت و احتمال مشخص و تحت درجه‌ای معین از آسیب‌پذیری‌های 
اقتصادی. اجتماعی و زیست محبطی تعریف می‌شود. ازاین‌رو ارزیابی و مدیریت ریسک خشکسالی به منظور کاهش 
اثرات حشکسالی و برنامه‌ریزی مناسب برای حشکسالی‌های آینده امری مهم است (وایت۲۰۱۲). از سوی دیگر 
تغیبرات اقلیمی؛ می‌تواند بر فراوانی و شدت وقوع واقعه خشکسالی تأثیرگذار باشد. براین اساس آگاهی از تغییرات 
وضعیت ریسک خشکسالی به دلیل تأثیر از تغییرات احتمالی اقلیمی. می‌تواند گامی مهم در کاهش خسارات ناشی از 
حشکسالی در بخش‌های مختلف (اقتصادیء اجتماعی» زیست‌محیطی و...) باشد. ریسک معمولاً به عنوان ترکیبی از 
احتمال رویداد یک واقعه (به‌عنوانمثال خشکسالی) و پیامدهای منفی آن شناخته می‌شود (شهید و بهروان» ۲۰۰۸ فن 
و همکاران" ۲۰۱۷). مقدار واکتش و پاسخی که محیط به مخاطره می‌دهد بستگی به‌شدت مخاطره و میزان ویژگی - 
های محیطی دارد. این ویژگی‌های محیطی را می‌توان به چهار دسته کلی عوامل اقتصادی» اجتماعی؛ فیزیکی و 
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زیست‌محیطی تقسیم نمود (برنامه استراتژی‌های بین‌المللی برای کاهش بلایا (151018)؛ ۰۲۰۰۷ ازاین‌رو تنها درک و 
پیش‌بینی مخاطره مهم نیست. بلکه عواقب احتمالی مخاطره. که بستگی به آسیب‌پذیری‌های محیطی و اجتماعی دارد 
نیز اهمیت دارند. اکثر سیستم‌های هشدار دهنده خشکسالی پایش و يا پیش‌بینی مخاطره را ارائه می‌دهند. اما اطلاعاتی 
در خصوص وجود و یا عدم وجود پیامدهای منفی خشکسالی ارائه نمی‌دهند (بچمایر و همکاران, ۲۰۱۷). اساسأه دو 
مرحله مهم در محاسبه مخاطره خشکسالی وجود دارد. اولین راه شناسایی وقایع حشکسالی با استفاده از شاخص‌های 
خشکسالی است» سپس احتمال خشکسالی محاسبه می‌شود (یوآن" ۲۰۱۳). برای بهبود مدیریت ریسک خشکسالی» 
شناسایی علل آسیب‌پذیری خشکسالی در یک منطقه حاص. گام مهمی است. پژوهش‌های مختلفی در ارتباط با 
تعیین ریسک خشکسالی چه در داخل ایران و چه در خارج از ایران انجام شده است. در عمده این مطالعات مخاطره 
خشکسالی براساس شاخص‌های خشکسالی و احتمال وقوع آن‌ها محاسبه شده است و برای محاسبه آسیب‌پذیری 
ریسک و يا حتی در بعضی موارد محاسبه ریسک» برخی از مطالعات از روش‌های سنجش از دور استفاده کرده‌اند 
(چوپر؛ ۲۰۰7 لین و چن, ۲۰۱۱ لوبیس و همکاران ۲۰۱۷) و در برخی دیگر از مطالعات از روش‌های آماری 
استفاده شده است. این مطالعات معمول عوامل اقتصادی و اجتماعی را به‌عنوان عوامل آسیب‌پذیر ذکر می‌کنند 
(تاسکاریس " ۲۰۰۷ پاندی و همکاران؛ ۲۰۱۰ کیم و همکاران؟ ۲۰۱۳ فن و همکاران"» ۲۰۱۷ پی و همکاران» 
۷ وو وهمکاران/ ۲۰۱۸ اپورو و همکاران" ۲۰۱۹ و دای و همکاران"؛ ۲۰۲۰). در هر دو دسته از مطالعات 
فوق برای بررسی میزان آسیب‌پذیری خشکسالی؛ استفاده از پرسش‌نامه و تکمیل آن توسط کشاورزان و مردم منطقه 
انجام نشده‌است. در حالی که در اين مطالعه, استفاده از پرسش‌نامه روش اصلی در محاسبه آسیب‌پذیری می‌باشد. 
محاسبه ریسک ابزاری مفید برای دید دقیق و مناسب از وضعیت خشکسالی در هر منطقه‌ای است. در واقع 
این دید. باعث کمک به مدیران برای اتخاذ تصمیم‌های بهتر برای سازگاری با خشکسالی در دوره‌های 
زمانی آتی می‌باشد. از این رو در این مطالعه سعی شده است برای پایش وضعیت ریسک خشکسالی زیر 
حوضه آبریز افین از داده‌های مشاهداتی و برای پیش‌بینی ریسک خشکسالی, از داده‌های ریزمقیاس شده 
اقلیمی حاصل از پروژه 0۳ استفاده گردد. در انتهاء با مدل‌سازی ریسک خشکسالی: میزان ریسک 
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۱۳۰ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


برای دوره‌های آتی محاسبه و با مقدار ریسک بدست آمده از داده‌های شبیه‌سازی شده اقلیمی مقایسه 


خواهد شد. 


۲- مواد و روش‌ها 

۱-۲ منطقه مورد مطالعه 

زیرحوضه افین دارای وسعت 1۵۵ کیلومترمربع و بخشی از حوضه پترگان- خواف می‌باشد که زیر 
مجموعه حوضه آبریز مرزی شرق ایران است. مساحت این حوضه ۳۲۹۸۰۰۰ هکتار است. این منطقه واقع 
در بین عرض‌های ۳۳ درجه و ۶۲ دقيقه شمالی تا ۲۳ درجه و ۵۸ دقيقه شمالی و ۵۲ درجه و ۲۵ دقیقه 
شرقی نا ۵۳ درجه و ۳۱ دقیقه شمالی می‌باشد. رودخانه دائمی ذر این زیر حوضه» رودخانه خوشبخت 
است. متوسط بارندگی در زیرحوضه افین حدود ۱۵۱/۶ میلی‌متر و متوسط دمای منطقه ۱۶/۷ درجه 
سانتی گراد می‌باشد (شکل ۱). اقلیم این زیرحوضه براساس اقلیم بندی دومارتن خشک و براساس 
اقلیم‌بندی آمبرژه خشک سرد می‌باشد. وجود اقلیم خشک و سرد زیرحوضه را مناسب برای کشت زرشک 
کرده‌است. اقتصاد روستای زیر حوضه افین از سالیان دور بر مبنای کشاورزی و محصول شاخحص زرشک 
بوده است. بیش‌ترین زرشک تولیدی ایران در اين منطقه به‌ئمر می‌رسد (سازمان جهاد کشاورزی شهرستان 


زیرکوه خراسان جنوبی۱۳۹۷). 


شکل ۱- موقعیت جغرافیایی زیرحوضه آبریز افین در کشور و استان خراسان جنوبی 


۲-۲ روش مطالعه 

این مطالعه با هدف بررسی و مدل‌سازی ریسک خشکسالی هواشناسی و هیدورلوژیکی و به صورت 
یک مطالعه موردی انجام شده است. برای محاسبه مخاطره ریسک خشکسالی‌های هواشناسی و 
هیدرولوژیک. در گام نخست از شاخص‌های خشکسالی 5۳۳۲1 و 8۳1 بعنوان شاخحص‌های خشکسالی 
هواشناسی و شاحص 591 بعنوان شاحص هیدرولوژیک استفاده شده است. دوره پایه برای این مطالعه از 
سال (۱۹۸۳-۲۰۱۵) می‌باشد. داده‌های مورد استفاده نیز شامل داده‌های ماهانه متوسط دمای حداکش 
متوسط دمای حداقل. متوسط دما در نزدیکی سطح زمین. بارش و داده‌های جریان رودخانه (رودخانه 
خوشبخت در زیرحوضه افین) می‌باشد. مقیاس‌های زمانی مورد استفاده. مقیاس فصلی زمستان (ژانویه- 
مارس)؛ بهار (آپریل - ژوئن) و پاییز (اکتبر- دسامبر» مقیاس ٩‏ ماهه ٩(‏ ماهه منتهی به ماه نوامبر فصل 
زراعی داشت و برداشت محصول زرشک) و ۱۲ ماهه منتهی به سپتامبر (سالانه) می‌باشد. 

از شانعص خشکسالی ٩۳1‏ به دلیل بررسی تأثیر همزمان چند پارامتر دمایی و بارش بعنوان پارامترهای 
موثر در خشکسالی. استفاده‌شد. برای محاسبه شاخحص خشکسالی 5۳۳1 نیاز به داشتن مقدار تبخیر- تعرق 
می‌باشد برای این منظور از روش تبخیر- تعرق هارگریوز- سامانی برای محاسبه تبخیر تعرق پتانسیل 
استفاده‌شد (سامانی ۵ ۲۰۰۰). 

شکاشن و همکاران" در سال ۲۰۱۷ شاخص خشکسالی 81 را ارائه نمودند. در واقع این شاخحص 
اصلاح‌شده‌ی شاخحص 8۳1 می‌باشد. در این روش بارش موثر در محاسبات شاخحص استفاده می‌شود. جهعت 
برآورد میزان بارش موثر در این تحقیق از روش تجربی 11518۳ استفاده شده است. این روش مناسب برای 
مناطق خشک و نیمه خشک می‌باشد. 
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۱۳۲ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


حدول ۱- کلاس‌های متفاوت شاخص خشکسالی 101.۳۲۲ و 5181 
(۵۵۵,2006 ۷۷ ک ۳۱12و :کلف ,۲6۵1:2017 11۵125 :]ناه ,6021,2010 مصهتمو-]۷۱۵۵۵ :۳۳ه) 


طبقات شاخص توصیف وضعیت 
بیشتر از ۲ ترسالی بسیار شدید 
۵ تا ۱/۹۹ ترسالی شدید 
۱ تا ۱/۹ ترسالی متوسط 
۹ -۰/۹۹ نرمال 
۱- تا -۱/٩‏ خحشکسالی متوسط 
۵- ۷ ۱/۹۹- خحشکسالی شدید 
کمتر از ۲- خشکسالی بسیار شدید 


شاخعص خشکسالی (1:06 ۴مم؟۳ )٩12702701760‏ همانند روش محاسبه شاخحص خشکسالی ۳7 
می‌باشد. با این تفاوت که در شاعص ]1 از داده‌های جریان رودخانه‌ای استفاده می‌شود )0 م12 ۱/۵4 
20017 

ماکان 

مخاطره خشکسالی به‌صورت احتمال وقوع خشکسالی در سطوح مختلف شدت خشکسالی تعریف 
شده است. برای محاسبه احتمال وقوع خشکسالی از رابطه ۳ استفاده شده است. 

باتوجه به اينکه سه شاخحص خشکسالی استفاده‌شده در این مطالعه دارای توزیع نرمال و همچنین 
طبقه‌بندی این خشکسالی‌ها مطابق جدول شماره ۱ است. بنابراین با استفاده از جدول 2 نرمال. احتمال 
وقوع هر طبقه خشکسالی محاسبه شد. مقادیر احتمال برای هر یک از طبقات مطابق جدول ۲ می‌باشد. 


حدول ۲- شدت خشکسالی. احتمال وقوع و وزن هر طبقه از خشکسالی برای شاخص‌های خشکسالی مورد 


بررسی 
وضعیت شدت هر طبقّه | احتمال وقوع هر طبقه | وزن هر طبقه 
خشکسالی خفیف ۹ تا ۰ ۳۱۹ ۱ 
خشکسالی متوسط ۹- تا ۱- ۰/۹۰۹ ۲ 
خشکسالی شدید ۹- تا ۱/۵- ۰/۹۵ ۳ 
خشکسالی خیلی شدید کمتر از ۲- ۳۲۳۹ ‌ 


۲۲ < ۲۷ ۶ ۳۵( 0۶۷ ۳۱ 


در رابطه ( ,1255 (۳0۵ احتمال وقوع هر طبقه کلاس خشکسالی و انآا شدت خشکسالی همان 
خشکسالی و ۷۷ وزن طبقه مورد نظر می‌باشد. 

۲-۲-۲ آسیب‌پذیری خشکسالی 

آسیب‌پذیری مجموعه‌ای از شرایط اقتصادی, اجتماعی و زیست‌محیطی است که حساسیت یک محیط 
به مخاطره را تعیین می‌کند (برنامه استراتژی‌های بین‌المللی برای کاهش بلایا (۲۰۰۷)15). برهمین اساس به- 
منظور شناسایی میزان تأثیر هریک از عوامل آسیب‌پذیری» پرسش‌نامه‌ای متشکل از ۳۵ سوال در سه بخش 
اقتصادی. اجتماعی و زیست‌محیطی طراحی شد. پرسش‌نامه در منطقه مورد مطالعه و به‌صورت تصادفی 
توسط کشاورزان تکمیل شد. تعداد نمونه لازم برای تکمیل پرسش‌نامه براساس روش‌های آماری و با 
استفاده از نرم افزار ۳۸۹5 ۵1055 4٩‏ عضو تعبین شد (هینتز؛ ۲۰۰۱). 

پس از تکمیل پرسش‌نامه برای تعیین مقدار آسیب‌پذیری منطقه. از روش‌های تحلیل عاملی اکتشافی " و 
تحلیل عاملی تاییدی* استفاده‌شد. درجه آسیب‌پذیری» عددی بین صفر تا یک بدست خواهد آمد. هرچه 
این مقدار به ۱ نزدیک‌تر باشد. به‌معنی بیشتر بودن میزان درجه آسیب‌پذیری محیط از حشکسالی می‌باشد. 

۳-۷-۲ نک عشکسالی 

پس از محاسبه مخاطره خشکسالی و آسیب‌پذیری آن, باتوجه به تعریف ریسک» می‌توان مقدار ریسک 


۱ / ۶۱ 


2-۲-۲ داده‌های اقلیمی آینده 

برای دستیابی به داده‌های هواشناسی مورد نیاز در دوره‌های آينده با توجه به میزان دسترسی به 
داده‌های پروژه 00۳ (دسترسی به خروجی‌های مدل‌های اقلیمی پروژه 00۳1۳ در سال ۱۳۹۱) و 
یز براساین آیتکه نایک نام پارامع‌های هزاشناسی وود نان دادههتای تاریش وک رح 
داشته‌باشند بنابراین تنها ۳ مدل (۱6۷۲2۲ سا 0۱-0۸ و ۲ هس10۲7۳0۳۲) انتخاب شد 
(جدول ۳). علت انتخاب این دو سناریو بررسی ریسک آینده در خوش‌بینانه‌ترین حالت و بدبینانه‌ترین 


۵00 تمافهد1۱1 ۲۵۲ رعماهتنه لقممتقصهام [ 
6 2 
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۱۳ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


دوره‌ی زمانی آینده در این تحقیق. سال‌های ۲۰۲۰-۲۱۰۰ می‌باشد که به سه دسته زمانی» آینده نزدیک (۲۰۶- 
20۳۰ انتیه میانی (۲۰۶۷-۲۰۱۷/۳) و آینده دور (۲۰۷-۲۱۰۰) تقسیم‌بندی سستل 6 استت: دلیل این تقسیم‌بندی» 
به‌منظور مقایسه‌ی مناسب تغییرات میزان ریسک خشکسالی میان دوره پایه (سال‌های ۱۹۸۹-۲۰۱۵) و دوه‌های 


زمانی آنی» می‌باشد. 


حدول ۳- مشخصات مدل‌های اقلیمی استفاده شده در این مطالعه 


تفکیک (رزولوشن) 7 
نام مدل شناسه موسسه کشور مرکز مدلسازی (یا گروه مدلسازی) 
افقی 
6 ۱۱2000۵1 امن تال ۷۲۵0۵16[ 
و 9۱ 9۱ ۳ ام - معنامنع۱۷۱6۲6۵۲0]0 ۲۹۵06۵۲006 
2۷5 9 ۵ 186 60 1560120۳06 06 ط۳۲۱۲00۵6۵ 
نوتمه ملق جه و۷۵66 
ان سارللن ۱9 شون هرا متا آممزوزطممع0) م۲۸ 
۸( نان ۸ ۱:9 
۶ ۶ ۰/1 ۳7 | ۳۳۳۸۵۲۲۲۲ | 1۳۵1220 صتا0عج۵0 ۳-۰۲۱۸۵۲۲۳۲۱ 
پس از استخراج داده‌ه ضریب همبستگی. ضریب تببین و ضریب نتش- ساتکلیف (مک کویین ! 


۹ به منظور بررسی توانایی مدل‌ها در شبیه‌سازی برای اه محاسبه شاد در مطالعات تغییر آب و 
هواء منابع مختلفی از عدم قطعیت وجود دارد که بر نتایج نهایی تأثیر می‌گذارد. یکی از منابع عدم قطعیت؛ 
خحطای مدل‌هاست است. به‌منظور از بین بردن این عدم قطعیت. یکی از روش‌ها ترکیب کردن مدل‌های 
اقلیمی است. از روش‌های رایج برای این منظور روش وزن‌دهی میانگین‌های مشاهداتی دما و بارش در 
قالب رویکرد بیسزین می‌باشد (معادله 1-۱ و 1-۲) (لشو و همکاران " ۷ هاشمی-آنا و همکاران " 


۵ طرحواره این روش برگرفته از مطالعه هاشمی - آنا و همکاران در شکل شماره ۲ مشاهده می‌شود. 
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در معادلات فوق. ‏ معرف بارش و ۲ معرف دما می‌باشد. ۸۳ اختلاف میانگین بارش شبیه‌سازی شده 
در دوره آینده از میانگین بارش مشاهداتی مقیاس‌های زمانی دوره پایه است (. * نیز تعداد مدل‌هامی- 


باش: 
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شکل ۲- طرحواره گام‌های روش بیزین براق آنالیز عدم قطعیت در این تحقیق (هاشمی - انا و اهمکاران) 


۲ ۲10۲ع)و۳۵ 
نامز( 
(م1 ,۳0۲۱۳ 


مره صفند۵۲٩۳‏ 


۵-۲-۲ شبیه‌سازی رواناب دوره‌های آتی 

برای دستیابی به داده‌های رواناب دوره‌های آتی از مدل 1۳۸61875 استفاده شده است. این مدل» یک 
مدل پارامتری بارش- رواناب است که در بیشتر مناطق آب و هوایی کاربرد دارد. باتوجه به داده‌های 
ورودی اندک و نیز صرف زمان کم به آسانی در بسیاری از حوضه‌های آبخیز این مدل به کار برده می‌شود 
(لیتوود و همکاران» 1۹۹۷). به‌منظور بررسی صحت‌سنجی نتایج از ضریب همبستگی و آزمون پارامتری 
احمل‌تتاه.] استفاده شد. 

۲-۲-۲ پیش‌بینی و مدل‌سازی ریسک خشکسالی 

به‌منظور پیش‌بینی ریسک خشکسالی از دو روش استفاده شد. در روش نخست (روش مستقیم). همانند 
دوره پایه. برای دوره‌های آتیء مقادیر خشکسالی» مقدار مخاطره و در انتها ریسک خشکسالی محاسبه 
می‌شود. در روش دوم. با کمک دو روش آماری» مدل‌سازی ریسک خشکسالی براساس ریسک خشکسالی 
دوره پایه» انجامش و سپس به کمک مدل‌های بدست آمده, ریسک خشکسالی برای سال‌های آتی؛ 

در مدل‌سازی ریسک خشکسالی باتوجه به اینکه مقدار درجه آسیب‌پذیری بدست آمده از پرسش‌نامه» 
در طی زمان‌های مختلف ثابت در نظر گرفته شده است؛ از این‌رو مدل‌سازی مخاطره انجام شده است. 

باتوجه به‌اینکه مخاطره خشکسالی براساس شاخص‌های خشکسالی بدست می‌آید و هر شاخحص 


1 116۷7000 6 1 


۱۳۹ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


ابتدا پارامترهای موثر در هر مخاطره خشکسالی (هواشناسی وهیدرولوژیک) شناسایی و سپس مدل مناسب 
برای پیش‌بینی بدست آمد. پس از پیش‌بینی مخاطره. با ضرب کردن مقدار آسیب‌پذیری که در قسمت قبل 
توضیح داده‌شد. مقدار ریسک خشکسالی برای سال‌های آتی محاسبه می‌شود. به منظور مدل‌سازی ریسک 
خشکسالی از دو روش مدل‌سازی 687655100 عاتامم0 300502260 (با تکرار ۰۵۰۰ ع1تاعصهتهم‌جمل ووعم1 
0 استفاده شد. 

«مندمی‌نوع3 علننمدنم 60مهتادا130: مدل‌های ر گرسیون معمولی. برای میانگین شرطی برازش داده می‌شود 
اما رگرسیون چندکی. برای چندک‌های شرطی بکار می‌رود. مهم‌ترین کاربرد رگرسیون چندکی شناسایی 
شکل توزیع متغیر وابسته‌ی مدل در سطوح گوناگون می‌باشد. ویژگی دیگر این روش مدل‌سازی, بوت 
استرپ کردن است. در واقع با این روش, تعداد نمونه‌ها را افزایش داده تا از اين طرق با ایجاد نمونه‌های 
فراوان» شرایط نمونه را به شرایط جامعه نزدیک کرده و با در نظر گرفتن تمامی حالات تشکیل نمونه از 
صحت برآورد ضرائب و برآورد فاصله اطمینان برای ضرایب. اطمینان بیشتری پیدا کرد (افرون؛ ۲۰۰۰). 

چندک‌های مورد بررسی در این مطالعه. چندک‌های ۰۲۵ ۵۰ ۷۵ و ٩٩‏ می‌باشد. 

ممنومع‌تعم1 متتاعجصمتهممع۱۱ وفع‌ما: ووعما یک مدل تعمیم یافته و متداول از رگرسیون محلی ابر کروتبوق 
چند جمله‌ای است که به رگرسیون میانگین متحرک نیز مشهور است.1۳۹۹ یک روش رگرسیون غیر 
پارامتری است. در رگرسیون پارامتری. شکل مشخصی برای پارامترهای مدل. در نظر گرفته می‌شود در 
حالی که در رگرسیون ناپارمتری شکل مشخصی برای پارامترهای مدل در نظر گرفته نمی‌شود. در روش 
ناپارمتری داده‌ها خود از میان مدل‌ها. بهترین مدل را برمی‌گزینند و این حالت باعث افزايش کارایی روش 
رگرسیون ناپارمتری نسبت به رگرسیون پارامتری می‌شود. 

از میان مدل‌های بدست آمده از چهار چندک استفاده شده (99 ,75 ,50 ,25-) در این تحقیق. چندکی که 
مدل مناسب‌تری براساس ضریب توافق ویلموت ((0 315۳۰ و ۷۸۲ را داشته‌است. انتخاب گردید. 


برای تعیین توانایی و مناسب بودن مدل‌ها از آماره‌های ذیل استفاده شده است. 
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ور تاش روا رش ون شاف مق ی له ترس بلاق قادیر سا هه دوس با شم بشید 
و یم 

براساس این روابط ریاضی توانایی و مناسب بودن مدل‌ها مشخص می‌شود. براساس آماره‌های عددی 
صمل طعملا ۳/55 و 2155۳ بهترین وضعیت. مقدار صفر (۰) و بدترین 0 است. براساس شاخص 
توافق ویلموت. مقدار یک (۱) بهترین وضعیت مقدار ۰ بدترین وضعیت است. 

۳ نتایج و بحث 

۱-۳- نتایج بررسی توانایی مدل‌های اقلیمی مورد بررسی و توسعه یک مدل ترکیبی 

برای بررسی توانایی پیش‌بینی مدل‌های اقلیمی استفاده شده. از ضریب همبستگی ((18 ضریب تبین 0۳2 
و ضریب نتش- سانکلیف (۵15 استفاده‌شد. نتایج این ضرایب حاکی از توانایی کامل مدل‌های دمایی 
(دمای حداقل, دمای حداکثر و متوسط دما در ارتفاع ۲ متری از سطح زمین) می‌باشد؛ اما هر سه مدل 
بارش استفاده شده در این مطالعهء نیازمند تصحیح می‌باشند» از این‌رو با استفاده از روش تغییر فاکتوره 
سعی در تصحیح مقادیر بارش برای هر سه مدل شد. پس از تصحیح سپس با استفاده از روش میانگین 
وزنی با رویکرد بیزین» مدل ترکیبی برای بارش و پارامترهای دمایی مشخص گردید. در واقع براساس این 
روش وزن هر یک از مدل‌ها تعیین و با استفاده از میانگین وزنی و بر اساس هر سه مدل, مقدار نهایی 
پارامتر مورد نظر محاسبه می‌شود. براساس نتایج این بخش. برای داده‌های بارش باتوجه به سناریوی 
5 مدل 11180-۲۸۲۲۲ بیش‌ترین ضریب را در مدل ترکیبی داشته است. براساس سناریو 30۳8.5 
مدل مناسب. در سه دوره آتی متفاوت است. اما به‌صورت کلی می‌توان گفت در آینده نزدیک مدل 81( 


در آینده میانی هر دو مدل 01501 و 1011580-۳۸17 به صورت یکسان در مدل تأثیر داشته‌اند و در آینده 


۱۳/۸ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


دور مدل 1013۳50-۳۸111 دارای بیش‌ترین تأثیر می‌باشد. برای هر سه پارامتر دمایی؛ و براساس هر دو 
سناریوی انتشار خوش‌بینانه و بدبینانه. مدل 0۳۲۲-۳۷۲2۷ بیش‌ترین وزن را در مدل ترکیبی داشته است. 

۲-۳- نتایج مدل 1۳1۸01815 

نتایج این بخش نشان داد که براساس دو سناریوی انتشار. مقدار رواناب در سال‌های آتی تغییر قابل 
ملاحظه‌ای نخواهد داشت. افزایش دما و بدنبال آن, افزایش احتمال وقوع بارش‌های حدی در منطقه مورد 
مطالعه شاید دلیلی برای تغییرات کم در مقدار رواناب در سال‌های آتی باشد. که افزایش دما در سال‌های 
آینده براساس خروجی مدل‌های 0018۳۳6 استفاده شده در این تحقیق, تأیید می‌شود. از دیگر دلایل 
احتمالی عدم کاهش رواناب. می‌توان به افزایش ذوب برف زیر حوضه‌های مجاور در اثر افزايش دما 


دانست که به صورت جریان‌های زیر سطحی و یا سطحی وارد منطقه می‌شوند. 


جدول 6- متوسط مقادیر رواناب. بارش سالانه. دمای حداقل. دمای حداکثر و دما در ارتفاع ۲ متری از سطح 
زمین در دوره پایه و دوره‌های آتی برمبنای 130۳4.5 و 130۳8.5 
پارامتر دوره پایه 10 
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۰/۰۰ | ۰۸۰۷ | ۰۸۷ | 85 
۱۲/۸۲ ۱ ۱۳۸۸۵۲۰ ۱ ۱۳۵۳ ۱ 45 ۱۳۶ 


روثاب (03/9) 


بارش سالانه (۲0۲0) 
8۵5 | ۱۳۸۸۹ | ۱۳۹۶ | ۱۲/۸۸ 


٩۸۷۹ 1 1/۶۲ | ۸۸۷۲ | 045 ۷/6‏ 
دمای حداقل ۵9 
85 | ۸۹۰ | ۱۰/۲۵ | ۱۱/۸۶ 


۲۲ 45 | :۲۳۰ | ۲۵۸۰۷ ۱ ۲۱/۵۸ 
دمای حداکثر (/ 
ک68 | ۲۳۵۸ 1 ۲۵/۸۰۲ 1 ۲۰۸۷۱ 


۱۷/۲۰ | ۱2/۷۶ | ۱۵/۸۶ 4,5 


دما در ارتفاع ۲ متری سطح زمین 0 ۱1/۸۷ 
کوصو ادا ۱۵ ۱۹۲۸ 


۳-۳- آسیب پذیری خشکسالی 
این عدد نشان می‌دهد که در حدود نیمی از ایلع زیرحوضه به دلیل عوامل اقتصادی. اجتماعی و 


می‌باشد. به صورت مجزا درجه آسیب پذیری برای هریک از عوامل اقتصادی. زیست محیطی و اجتماعی به 
ترتیب عبارتند از ۰/۳۰ ۰/۳۶ و صفر. براین اساس. عامل زیست‌محیطی سهم بیشتری در آسیب‌پذیری 
زیرحوضه افین را دارا می‌باشد. در مقابل عامل اجتماعی تاثیری در آسیب‌پذیری منطقه ندارد. سیب پذیری 
برای دوره پایه ۰/۵۳ و برای دوره آتیء بازه ۰/۳۰ تا ۰/۷۰ در نظر گرفته شد. 
(جمتاجز06 ۲0ح0صمای + بوانانطدصهان۷) 

2-۳- ریسک خشکسالی براساس دو روش مستقیم و استفاده از مدل‌های پیش‌بینی 

نتایج این بخش نشان می‌دهد که براساس شاخحص خشکسالی 081 در دوره‌های آتی شدت وقایع 
ریسک خشکسالی نسبت به دوره پایه نسبتاً بیشتر می‌شود بدین ترتیب که تعداد و شدت وقایع خشکسالی 
(شامل چهار وضعیت خشکسالی خیلی شدید. شدید. متوسط و خفیف) در دوره آینده دور نسبت به دو 
دوره آینده میانی و نزدیک براساس سناریوهای انتشار بدبینانه و خوش‌بیبنانه بیشتر می‌شود. همچنین نتنایج 
وقوع خشکسالی‌ها و شدت‌های آن‌ها در دوره‌های آتی نشان داد که شدت خشکسالی‌ها متأثر از سناریوی 
انتشار نیز می‌باشد بدین ترتیب که با افزایش میزان انتشار گازهای گلخانه‌ای. شدت و یا حتی تعداد وقایع 
خشکسالی‌ها افزایش می‌یابد که این مطلب در نتایج بدست آمده از شاحص خشکس‌الی 5۳۳۲ و ۰80 کامل 
نمایان است» چرا که در این شاخص علاوه بر بارش, پارامترهای دمایی نیز دخیل هستند. با افزايش گازهای 
گلخانه‌ای حبس شدن گرمای (مادون قرمز) منعکس شده از سطح زمین افزایش و درنتیجه باعث بالارفتن 
دماهای سطحی خواهد شد. 

نتایج ریسک خشکسالی نشان می‌دهد که در فصل زمستان بیشترین تغییرات تعداد دفعات خشکسالی 
نسبت به دوره پایه حادث خواهد شد. براساس شاخص ٩81‏ نیز ریسک‌های خیلی شدید ناشی از 
خشکسالی‌های خیلی شدید در هیچ یک از مقیاس‌های زمانی مورد بررسی در سه دوره پراساس هردو 
سناریوی انتشار. مشاهده نشد. همچنین مشخص شد که در دوره آینده نزدیک» فراوانی وقوع ریسک 
خشکسالی‌های شدید و متوسط نسبت به دوره پایه و دو دوره زمانی آینده میانی و دور بیشتر خواهد بود. 

باتوجه به شاحص ٩۳۲1‏ و ]5 که چهار پارامتر متوسط دما متوسط دمای حداقل. متوسط دمای 
حداکثر و بارش برای محاسبه‌ی مقدار آن لازم می‌باشد براساس روش‌های مدل‌سازی و انتخاب بهترین 
مدل پارامترهای موْثر در برآورد مخاطره خشکسالی در هریک از مقیاس‌های زمانی مورد مطالعه مشخص 


و پیش‌بینی میزان مخاطره و پس از آن, محاسبه وت کنخ صورت پذیرفت. 


۱۳۰ جغرافیا و مخاطرات محیطی شمارة سی و سوم 


به‌منظور صحت سنجی مدل‌هاء ۱۰ سال آخر دوره پایه برای این منظور استفاده شد. بدین ترتیب که ۲۳ 
سال نخست (۱۹۸۳-۲۰۰۵) برای انتخاب مدل و ۱۰ سال پس از آن (۲۰۰۳-۲۰۱۵) برای صحت سنجی 
لها انتفاده شک 

باتوجه به نتایم بدست آمده. همراستایی قوی بین مقادیر مخاطره حاصل از مدل‌های پیش‌بینی و مقادیر 
مخاطره مشاهده شده براساس سه شاخحص خشکسالی مورد مطالعه و براساس هر دو سناریو دیده می‌شود. 

پراساس روش 808 و با توجه به آماره‌های شاحص توافق ویلموت ((0 1۷5۳۰ و ۰۷۸۲ مناسب‌ترین 
مدل در میان چهار مدل (99, 25,50,75-). مدل با چندک ۵۰ (چارک میانه) در بیشتر مقیاس‌هاو در دو 
شاخص خشکسالی هواشناسی وبرای شاخص خشکسالی هیدرولوژیک. مدل چندک ۸٩‏ تعیین شد (در 
جداول ۵ * و ۷ چارک مناسب تعیین شده برای هر مقیاس به تفکیک شاخص‌ها ذکر شده است). در اینن 
قسمت نیز ابتدا پارامترهای موثر بر شاحص‌های ٩۳۳1‏ و 08 مشخص و سپس مدل‌سازی انجام شد. 
پراساس شاخحص ٩۳۳۲‏ در ۶ چندک. بارش پارامتر موثر شناخته شد. تنها در مقیاس زمانی ٩‏ ماهه متوسط 
دمای حداقل و متوسط دما در نزدیکی سطح زمین موثربوده است. با توجه به شاخص 18 در فصل 
زمستان علاوه بر بارش متوسط دمای حداکثر متوسط دمای حداقل و متوسط دما در نزدیکی سطح زمین 
مثر بوده است. در سایر مقیاس‌ها. بارش. پارامتر مزثر برای پیش‌بینی می‌باشد. 

براساس روش رگرسیون ناپارمتری 10655 در دو شاخحص خشکسالی ٩۳۳1‏ و 8۳ با توجه به اینکه 
چهار متغیره می‌باشد. به‌منظور شناسایی پارامتر موثر در مدل. نخست از روش رگرسیون ناپارامتری ساده 
بوت استرپ شده (با تکرار 4۵۰۰ پارامتر موثر شناسایی و سپس مدل 10655 برای پارامتر موثر اجرا شد. 
در تمامی پنج مقیاس زمانی براساس دو شاحص خشکسالی 5۳۳۲ و 1 پارامتر بارش تنها پارامتر موثر 
شناخته شد. به منظور انتخاب بهترین مدل پیش‌بینی براساس دو روش مذکور از آماره‌های عددی تعیین 
به‌منظور تعبین توانایی و مناسب بودن مدل‌ها استفاده شد. نتایج این بخش در جداول ٩۰۵‏ و ۷ گزارش 
شده است. همانطور که از جداول مشسخص است. در بیشتر مقیاس‌های زمانی و برای سه شاخحص 
خشکسالی ٩۳۳۲‏ . اه و 5۳1 مدل رگرسیون ناپارامتری 10665 بیش‌ترین توانایی را در پیش‌بینی مخاطره 
داشته است. 

براساس نتایج این بخش نیز می‌توان گفت که روش رگرسیون ناپارامتری 10085 در شاخص ٩1‏ اگرجه 
در بیشتر مقیاس‌ها بعنوان مدل مناسب شناسایی نشده است. اما مقدار برآورد میزان مخاطره برای سال‌های 


آتی را مانند روش مستقیم. با دامنه بیشتری از تغییرات نشان می‌دهد. 


حدول ۵- شاخص توافق ویلموت (۰)0 ۸5 ۷۸0 ۷۲5۲۲ و :1151 برای مدل‌های ریسک شاخص 5۳۲۲ 


گام زمانی مدل | ۱۲6۲۴ | ۲۸0 | لها ۰ 4 


094 | 3 003 0 004 10۳55 
072 ۱ 61 0.06 1 01 0.08 ۰. 30۵-50 
078 ۱ 0.04 1 1 007 10355 
085 | 9 005 0 0.06 ۰. 30۵-50 
095 | 3 003 0 004 10055 
0۳21۹۲ ۰ ۰.4 0.07 1 01 0.08 ۰. 30۵-50 


اماه منتهی به فصل برداشت | 10۳۴855 009 001 009 ۵ | 052 


زرشک وفو5080 .012 | 001 .| 01 .| 238 [ 062 
۳3 5 003 0 | 003 | 034 | 096 
(25-)08ظ 02 0 1 004 | 048 | 0.89 


۵-۳- نتایج مقایسه میان مقادیر ریسک خشکسالی پیش‌بینی شده توسط مدل‌های پیش‌بینی با رسک خشکسالی از 
روش مستقیم در سال‌های آتی 

در این بخش به مقایسه میان نتایج ریسک خشکسالی محاسبه شده به روش مستقیم تا نات 
خشکسالی‌های پیش‌بینی شده توسط روش‌های مدل‌سازی پرداخته شده است. برای این منظور از مقدار 10 
۷20 آزنون 6۲ او ضریب میت کین استفاده شده است. (حداول ۸ و 4) 

باتوجه به نتایجی که در جحداول ۸و ٩‏ ارائه شده است. می‌توان دریافت که در محاسبه ریسک 
شفی‌کسای باق ال هبای انیس براساس روش مستقیم و روش مدلسازی با استفاده از مدل مناسب؛ 
همراستایی قوی مشاهده می‌ شود (باتوجه به بالا بودن مقادیر ضریب همبستگی). براساس دو سناریوی 
انتشار خوش بینانه و بدبینانه. در مقیاس زمانی سالانه و فصل بهار مقدار کمتر همراستایی مشاهده می‌ شود. 
براساس نتایج آزمون 1-000601 نیز می‌توان گفت که اختلاف معنی‌داری در مقدار برآورد ریسک خشکسالی 
براساس هر دو روش وجود ندارد. تنها در شاعص 7 در مقیاس‌های زمانی بهار و ٩‏ ماهه منتهی به 
مشاهده می‌شود که دلیل این امر می‌تواند به دلیل پراکندگی بیشتر بارش در این مقیاس‌ها و توام شدن آن با 

تر ان یی هون سای وه اه ایا تس هی استا ها تالا در ان اس ماف رن کت 
که مدل‌ها با اعمال اصلاحاتی مناسب هستند. 


گام زمانی 


ماه منتهی به فصل برداشت زرشک 


سالانه 


جغرافیا و مخاطرات محیطی 


مدل 


10055 
(80۳)2-75 
10۳55 
380100250 
10۳55 
13801)0-50 
10۳55 
(50-)0۲ظ۳ 
10۳055 
80100-250 


(5۳, 
0۳05 
011 
0۳05 


(۱9۳ 
)01 
)(01 
001 


001 
001 


۱/۳۹۱ 


0۳05 
009 
003 
006 
003 
0۳05 
002 
003 
0۳05 
006 


۱/۳۹۹۹ 


)1 
0۳87 
0۳69 
144 
)1 
129 
0۳24 
044 
0۳63 
00 


شمارة سی و سوم 


حدول 1- شاخص توافق ویلموت (۰0 1۸۲ ۸ ۷15۲۳ و 1۱15۲ برای مدل‌های ریسک شاخص ۰1۲ 


0 


012 
035 
0۳12 
0853 
093 
0۳12 
096 
092 
060 
062 


حدول ۷- شاخص توافق ویلموت (۰0 ۷/۸1۲ ۸۲۸ ۱151۳ و 13۷15 برای مدل‌های ریسک شاخص 51٩1‏ 


گام زمانی 
پاییز 


بهار 


ماه منتهی به فصل برداشت زرشک 


سالانه 


مدل 
۱992 
(۲01)<99 
۱992 
(0۵۲)0-099ظ۳ 
۱992 
۷6 ۱۰۵ 
لد ۱9 
۷6 (۱۰/۵ 


10۳055 
1801)9-99( 


(5۳ 
02 
012 
001 
011 
018 
0۳05 
019 
007 


0۳15 
0۳05 


(9۳ 
0۳04 
002 

0 
001 
003 
001 
003 
001 


002 
)(01 


۱/۳۹۱ 
018 
011 
001 

01 
012 
006 
015 
006 


011 
10 


۱/۳۵۹۹ 
109 
0۳66 
0۳04 
06 
0۳68 
033 
092 
040 


013 
04 


0 
017 
041 
099 
049 
051 
08 
040 
06 


05 
0۳617 


حدول ۸- ۲ مقدار آزمون 1[ و ضریب همبستگی 0 میان مقادیر ریسک خشکسالی آینده برمبنای روش 


مستقیم و مدل مناسب سه شاخص خشکسالی 5۳۲1 1801 و 5181 براساس 1304.5 


گام زمانی 


59۳۳ 
ِ مدل مناسب 
۱-06 

ِ 5[ 
ح 5[ 
و _ تسیا 
ت 4( 
هون نما 


(» 
ِ مدل مناسب 
۴-6( 
توس ما 
ك ۹( 
10695 
6 | ععما 
بو | ۵« 


ال 


4 
[5 
4 
04 


(4 


حدول 9- ۲ مقدار آزمون ۲[ و ضریب همبستگی 003 میان مقادیر ریسک خشکسالی آینده برمبنای روش 


مستقیم و مدل مناسب سه شاخص خشکسالی 5۳۲1 1801 و 5181 براساس 1308.5 


گام زمانی 9۳ 0 4 
مدل مناسب ِ مدل مناسب ِ مدل مناسب 
۱-06 ۱-06 ۱-6 
: 091 015 013 
زمستا 
زمستان 0 1695[ 0005 1695 0 04 
014 ۳ 094 
بای 
پاییز 0 وومم ][ 005 04 0 وومم [ 
بهار ِ 5و6 1 ِ 1695 اب 64 
001 024 0 
۱ ۱ : 035 015 0/6 
ماه مت اشت زرشک 
منتهی به فصل برداشت زر 0 04 0 10255 0 ۳0 
0۳68 0036 014 
سالا 
3 017 ی 0 ِ 


- نتیجه گیری 

این مطالعه, با هدف شناسایی ریسک خشکسالی هواشناسی و هیدرولوژیک و توسعه مدلی برای 
پیش‌بینی آن به‌منظور پیشگیری از بروز خسارات شدید و جدی به‌صورت موردی در زیر حوضه آبریز افین 
انجام شد. برای این‌منظور نخست از شاخص‌های خشکسالی هواشناسی و هیدرولوژیک بعنوان ابزاری 
برای اندازه‌گیری شدت خشکسالی‌ها و به دنبال آن محاسبه میزان مخاطره ریسک خشکسالی‌های استفاده 
شد. آسیب‌پذیری نیز با کمک پرسش‌نامه از کشاورزان منطقه تهیه و در انتها میزان ریسک خشکسالی برای 
دوره زمانی پایه محاسبه گردید. برای پیش‌نگری میزان ریسک خشکسالی برای دوره‌های زمانی آتی نیاز 
به داده‌های اقلیمی و هواشناسی دوره‌های زمانی آتی بود که از این‌رو از مقادیر خروجی پروژه اقلیمی 
۴ و نیز مدل بارش- رواناب 111۸018۳5 استفاده شد. نتایج بررسی سه مدل اقلیمی صحت و 
درستی سه پارامتر دمایی را تأیید نمود ولیکن پارامتر بارش نیز با توجه به اقلیم نیمه عشک و پراکندگی 
بارش در منطقه افین» پس از اعمال تصحیحات. قابل استفاده گردید. در پیش‌نگری‌ها از مقادیر رواناب 
حاصله از مدل 111۸08۳5 استفاده شد. چراکه در مقایسه میان پارامتر رواناب براساس خروجی مدل‌های 
اقلیمی 6080۳ و مدل بارش- رواناب 11۸08۳5 نتایج بدست آمده از روش 111۸018۳5 دقت بالاتری 
را نشان داد. رزاقیان و همکاران (۱۳۹۵) و پور محمدی و همکاران (۱۳۹7) کاهش رواناب را برای سال- 
های آتی با استفاده از دو روش مدل‌سازی 111۸08۳85 و الگوریتم احتساب خاک پیش‌بینی نموده‌اند. در 


۱۳ جغرافیا و مخاطرات محیطی شمار؛ سی و سوم 


سال‌های آتی پیش‌بینی نموده‌اند؛ اما در تحقیق حاضر نتایج پیش‌نگری مقدار جریان سطحی براساس مدل 
رواناب برای سال‌های آتی است. 

سال‌های آتی شد. براین اساس و به منظور افزایش دقت در پیش‌بینی میزان ریسک خشکسالی‌ها» از دو 
روش مدل‌سازی استفاده شد. به‌منظور صحت‌سنجی هر یک از روش‌های مورد استفاده, برآورد مخاطره 
برای ده شا اخیر (۲۰۰۹-۲۰۱۵) انجام شنك: نتایج بخش صحت‌سنجی مدل‌ها. نشان داد که:ر فرسیون 
پس از این مدل. مدل رگرسیون نیمه پارامتری چندکی بوت استرپ شده مناسب تشخیص داده شد. برای 
پارامتر رواناب» روش مدل‌سازی رگرسیون نیمه پارامتری بوت استرپ شده با چندک ٩‏ در بیشتر 
مقیاس‌های زمانی با دقت بالاتری توانسته است. برآورد میزان مخاطره ۳ داشته باشد. بر اساس نتایج این 
بخش می‌توان گفت که روش‌های مدل‌سازی امکان برآورد و پیش‌نگری میزان ریسک خشکسالی‌ها را 
نگری مقادیر مخاطره برای دوره‌های زمانی آتی با استفاده از مدل مناسب انتخاب شده در هر مقیاس زمانی 
و براساس هر شاخص خشکسالی. با مقادیر مخاطره دوره آتی که به روش مستقیم بدست آمد. همراستایی 
قوی مشاهده شد اما اختلاف میانگین میان این دو مقدار معنادار شده است. باتوجه به اينکه مقادیر مخاطره 
روش مستقيم خود براساس پیش‌نگری متغیرهای هواشناسی هستند و نتایج مدل‌ها نیز براساس پیش‌نگری- 
های متغیرهای هواشناسی و عدم قطعیت‌های ناشی از انتخاب مدل هستند. شاید دلیلی بر معنی‌دار شدن 
میانگین‌ها باشد. دوره آماری طولانی مورد مقایسه (دوره ۸۱ ساله) می‌تواند عامل دیگر ایجاد اين اختلاف‌ها 
باشد. دلیل اهمیت دقت در انتخاب مدل‌های مناسب برای پیش‌نگری ریسک به دو دلیل است؛ نخست 
خشکسالی بر محیط است که اين پدیده را تبدیل به یک خطر می‌کند. دوم؛ اگرچه نمی‌توان مانع از وقفوع 
حوادث طبیعی چون خشکسالی شد اما می‌توان با داشتن آگاهی از وضعیت ریسک خشکسالی برای 
سال‌های آتی. تمهیدات لازم برای سا زگاری با این ید یده را به‌منظور به حداقل رساندن اثرات منفی واقعه 
خشکسالی فراهم ساخت. 
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